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We present measurements of three-dimensional correlation functions of like-sign low transverse momentum 
kaon pairs from y / 5/^v=200 GeV Au+Au collisions. A Cartesian surface-spherical harmonic decomposition 
technique was used to extract the kaon source function. The latter was found to have a three-dimensional 
Gaussian shape and can be adequately reproduced by Therminator event generator simulations with resonance 
contributions taken into account. Compared to the pion one, the kaon source function is generally narrower and 
does not have the long tail along the pair transverse momentum direction. The kaon Gaussian radii display a 
monotonic decrease with increasing transverse mass mj over the interval of 0.55<mx<1.15 GeV/c . While the 
kaon radii are adequately described by the my-scaling in the outward and sideward directions, in the longitudinal 
direction the lowest mj value exceeds the expectations from a pure hydrodynamical model prediction. 
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Analysis of the data collected at the Relativistic Heavy Ion 
Collider (RHIC) has resulted in the discovery of strongly in- 
teracting, almost perfect fluid created in high energy nucleus- 
nucleus collisions JH-0]. Lattice calculations predict that 
the transition between normal nuclear matter and this new 
phase is a smooth cross-over fl. This is consistent with 
the absence of long source lifetimes which would indicate a 
first-order phase transition . Moreover, analysis of three- 
dimensional (3-D) two-pion correlation functions, exploiting 
the novel technique of Cartesian surface-spherical harmonic 
decomposition of Danielewicz and Pratt 13, HI], revealed sig- 
nificant non-Gaussian features in the pion source function 13] ■ 
Furthermore, the extraction of the shape of the pion source 
function in conjunction with model comparisons has permit- 
ted the decoupling of the spatio-temporal observable into its 
spatial and temporal aspects, and the latter into source lifetime 
and emission duration. However, an interpretation of pion cor- 
relations in terms of pure hydrodynamic evolution is compli- 
cated by the significant contributions of resonance decays. A 
purer probe of the fireball decay could be obtained with kaons 
which suffer less contribution from long lifetime resonances 
and have a smaller rescattering cross-section than pions. The 
lower yields, however, make it difficult to carry out a de- 
tailed 3-D source shape analysis of kaons. A one-dimensional 
kaon source image measurement was recently reported by the 
PHENIX Collaboration 1 1 Oil - This measurement, however, 
corresponds to a fairly broad range of the pair transverse mo- 
mentum 2kj which makes the interpretation more ambiguous. 
In particular, information about the transverse expansion of 
the system, contained in the kj dependence of the emission 
radii, is lost. The one-dimensional nature of the measurement 
has also less constraining power on model predictions than 
would be available from a 3-D measurement. 

This paper presents a full 3-D analysis of the correlation 
function of mid-rapidity, low transverse momentum like-sign 
kaon pairs. The technique used in this paper is similar to that 
employed in the first 3-D extraction of the pion source func- 
tion It involves the decomposition of the 3-D kaon corre- 
lation function into a basis of Cartesian surface-spherical har- 
monics to yield coefficients, also called moments, of the de- 
composition which are then fitted with a trial functional form 
for the 3-D source function. The latter is then compared to 
models in order to infer the dynamics behind the fireball ex- 
pansion. 

The presented data from Au+Au collisions at ^/snn = 200 
GeV were taken by the STAR Collaboration during the year- 
2004 and 2007 runs. A total of 4.6 million 0-20% central 
events were used from year 2004, and 16 million 0-20% cen- 
tral events from year 2007. We also analyzed 6.6 million 
0-30% central events from the year 2004 run to compare to 
the previously published PHENIX kaon measurements [10]. 
Charged tracks are detected in the STAR Time Projection 
Chamber (TPC) iflTll . surrounded by a solenoidal magnet pro- 
viding a nearly uniform magnetic field of 0.5 Tesla along the 
beam direction. The TPC is used both for the tracking of 



charged particles at midrapidity and particle identification by 
means of ionization energy loss. The z-position of the event 
vertex is constrained to be \z\ < 30 cm. 

The 3-D correlation function C(q) = Af sam e(q) /Admixed (q) 
is constructed as the ratio of the 3-D relative momentum 
distribution, A^ same (q), for K + K + and K K pairs in the 
same event to that from mixed events, -/V m j xe d(q)- Here, q = 
(Pi — P2)/2, where pi and p2 are the momentum 3-vectors 
of the particles in the Pair Center-of-Mass System (PCMS). 
The non-commutativity of the Lorentz transformations along 
non-collinear directions demands that the Lorentz transfor- 
mation from the laboratory frame to the PCMS is made by 
first transforming to the pair Longitudinally Co-Moving Sys- 
tem (LCMS) along the beam direction and then to the PCMS 
along the pair transverse momentum. C(q) is flat and normal- 
ized to unity over 60 < |q| < 100 MeV/c. 

To obtain the moments, the 3-D correlation function C(q), 
is expanded in a Cartesian harmonic basis HI 

C(q)-l=*(q)= £ <...« I («X...a I («q) . d) 
/,a[ ...a; 

where / = 0, 1,2,..., a, = x,y orz, and a (£2q) are Carte- 
sian harmonic basis elements (i2 q is the solid angle in q 
space). R l a a Aq), where q is the modulus of q, are Carte- 
sian correlation moments, 

<...«,(?) = / ^...^^(q). (2) 

The coordinate axes x-y-z form a right-handed out-side-long 
Cartesian coordinate system. They are oriented so that the 
z-axis is parallel to the beam direction and x points in the di- 
rection of the pair total transverse momentum. 

Correlation moments can be calculated from the measured 
3-D correlation function using Eq. Even moments with 
/ > 4 were found to be consistent with zero within statistical 
uncertainty. As expected from symmetry considerations, the 
same was also found for odd moments. Therefore in this anal- 
ysis, the sum in Eq. (1]) is truncated at I = 4 and expressed in 
terms of independent moments only. Up to order 4, there are 6 
independent moments: R°, R^ x , /??L R xxxx , Rim an d^ De- 
pendent moments are obtained from independent ones 0, |8|] . 

These independent moments were extracted as a function 
of q, by fitting the truncated series to the measured 3-D cor- 
relation function with the moments as free parameters of the 
fit. The statistical errors on the moments reflect the statistical 
error on the 3-D correlation function. In order to estimate the 
effect of systematic errors, the 3-D correlation function and 
associated moments were obtained under varying conditions 
including nominal vs. reverse magnetic field, year 2004 vs. 
year 2007 data, positively vs. negatively charged kaon pairs 
and varying kaon sample purities. Although the variations did 
not introduce any observable systematic deviation in the cor- 
relation moments, they have some effect on the parameters of 
the 3-D Gaussian fit of Eq. ©. 
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FIG. 1. Independent correlation moments R l (q) for orders / = 0,2,4 
for mid-rapidity, low transverse momentum kaon pairs from the 20% 
most central Au+Au collisions at ^/snn=200 GeV. Panel (a) also 
shows a comparison between R®(q) and R(q). The error bars are 
statistical. The solid curves represent results of the Gaussian fit. 
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Figure Q] shows the independent correlation moments 
' ai tt; up to order 1=4 (open circles) for mid-rapidity 
(|y|<0.5), low k T (0.2</t T <0.36 GeV/c) kaon pairs produced 
in the 20% most central Au+Au collisions at x /saw=200 GeV; 
kj is half the transverse momentum of the pair. In panel (a), 
R°(q) is shown along with the one-dimensional correlation 
function R(q) = C(q) — 1 (solid circles); both represent angle- 
averaged correlation functions, but R°(q) is obtained from the 
3-D correlation function via Eq. d2} while R(q) is evaluated 
directly from the one-dimensional correlation function. The 
data points have been corrected for the effect of track mo- 
mentum resolution. The agreement between R°(q) and R(q) 
attests to the reliability of the moment extraction technique. 
Figures |TJb)-(f) show that while second moments are already 
relatively small compared to their errors, fourth moments are 
insignificant without any visible trend. This further justifies 
truncating Eq. ([TJ at / = 4. 

The probability of emitting a pair of particles with a pair 
separation vector r in the PCMS is given by the 3-D source 
function S(r). It is related to the 3-D correlation function C(q) 
via a convolution integral as 



C(q)-l=tf(q) = / (|(Kq,r)| 2 -l)S(r)</r 



(3) 



where the relative wave function (j)(q,r) serves as a six- 
dimensional kernel, which in our case incorporates Coulomb 
interactions and Bose-Einstein symmetrization only iH. 



Strong final state interactions are assumed to be negligible due 
to the small s-wave scattering length (~0. 1 fm) of two iden- 
tical kaons lfl3ll . Hence, no correction to the measured corre- 
lation function for Coulomb and other final-state interaction 
effects is required. Analogously to Eq. ([T), the source func- 
tion can be expanded in Cartesian harmonics basis elements 
as 5(r) = 

,<xi...ot( ^ai...a; (''Mo 1 ...0| (^r)- Eq. can then be 
rewritten in terms of the independent moments Q [8| . 

The 3-D source function can be extracted by directly fitting 
the 3-D correlation function with a trial functional form for 
S(r). Since the 3-D correlation function has been decomposed 
into its independent moments, this corresponds to a simulta- 
neous fit of the six independent moments with the trial func- 
tional form. A 4-parameter fit to the independent moments 
with a 3-D Gaussian trial function, 

(4) 

yields a % 2 /ndf =1.7. The correlation strength parameter 
X represents the integral short-distance contribution to the 
source function (14]. Figure Q] shows the fit as solid curves, 
making it evident that the quality of the fit is predominantly 
driven by the relatively small errors of R°(q). The values of 
the Gaussian radii and the amplitude (R x ,R y ,R z , X) are listed 
in Table |D 

Figure |2ja)-(c) illustrates the kaon correlation function 
profiles (circles) in the x, y and z directions (C(q x ) = 
C(q x ,0,0), C(qy) = C(0,q y ,0) and C(q z ) = C(0,0,q z )) re- 
spectively, obtained by summation of the relevant correlation 

terms C l ai . Ml (q) = 8/, + #a 1 ...a i _(?Ma 1 ...a,( ii q) U P to order 
1=4. The peak at qm20 MeV/c is coming from an expected 
interplay of Coulomb repulsion at q— >0 and Bose-Einstein en- 
hancement. The correlation profiles from the data are well 
represented by the corresponding correlation profiles from the 
Gaussian fit (line). Hence, the trial Gaussian shape for the 
kaon source function seems to capture the essential compo- 
nents of the actual source function. 

Figure Oa)-(c) depicts the extracted source function pro- 
files in the x, y and z directions (S(r x ) = S(r x ,0,0), S(r y ) = 
S{0,r y ,0) and S(r z ) = S(0,0,r z )) obtained via the 3-D Gaus- 
sian fit (dots) to the correlation moments. The two solid 
curves around the Gaussian source function profiles represent 
the error band arising from the statistical and systematic errors 
on the 3-D Gaussian fit parameters, as well as the uncertainty 
from the source shape assumption estimated using a double- 
Gaussian fit. Note that the latter becomes important for large 
r values only. 

The source function profile S(r y ) in the side direction re- 
flects the mean transverse geometric size of the emission 
source, while the source lifetime determines the extent of the 
source function profile S(r z ) in the long direction. Being in 
the direction of the total pair transverse momentum (hence the 
direction of Lorentz boost from the LCMS to PCMS frame), 
the source function profile in the out direction S(r x ) is char- 
acterized by the kinematic Lorentz boost, mean transverse ge- 
ometric size as well as source lifetime and particle emission 
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FIG. 2. Kaon correlation function profiles (circles) for mid-rapidity, 
low transverse momentum kaon pairs from the 20% most central 
Au+Au collisions at ,/^viv=200 GeV (a) C(q x ) = C(q x ,0,0), (b) 
C(q y ) = C(0,q y ,0) and (c) C(q z ) = C(0,0,q z ) in the x, y and z di- 
rections. The curves denote the Gaussian fit profiles. 



duration. In order to disentangle these various contributions, 
the Monte Carlo event generator Therminator lll5l is used to 
simulate the source breakup and emission dynamics. 

The basic ingredients of the Therminator model employed 
in the analysis are: (1) Bjorken assumption of longitudinal 
boost invariance, (2) Blast-Wave (B W) expansion in the trans- 
verse direction with transverse velocity profile semilinear in 
transverse radius p 0: v r (p) = (p/p m ax)/(p/p max + Vt J> 
where v ( =0.445 is obtained from BW fits to particle spec- 
tra fnll . (3) after a proper lifetime x, a thermal emission of 
particles takes place from the source elements distributed in a 
cylinder of infinite longitudinal size and finite transverse di- 
mension p max . At the point of source breakup, all particle 
emission is collectively viewed as happening from a freeze- 
out hypersurface defined in the p-x plane as x = To +ap. 
Hence, particles which are emitted from a generic source el- 
ement with coordinates (z,p) will have emission time t in the 
laboratory frame given by f 2 = (xo + ap) 2 + z 2 ■ 

Note that the BW mode of fireball expansion means that 
a = III 80 making x independent of p. Each source element is 
thus defined by only one value of the proper breakup time X = 
Xo and all particle emission from this source element happens 
instantaneously in the rest frame of the source element and 
the proper emission duration Ax is set to 0. Later on, we will 
also discuss another choice for parameter a which was used to 
describe the pion data 131 . 

Using a set of thermodynamic parameters previously tuned 



FIG. 3. Kaon source function profiles extracted from the data (solid 
circles with error band) and 3-D pion source function (squares) from 
PHENIX together with Therminator model calculation for kaons 
with indicated parameter values (triangles). 

to fit charged pion and kaon spectra lfl8ll . mid-rapidity kaon 
pairs at low &t were obtained from Therminator with all 
known resonance decay processes on and off. They were then 
boosted to the PCMS to obtain source function profiles for 
comparison with corresponding profiles from the data. 

Figure |3ja)-(c) indicates that the 3-D source function gen- 
erated by the Therminator model in the BW mode (solid tri- 
angles) with Xo=8.0±0.5 fm/c, p max =9.0±0.5 fm and other 
previously tuned parameters lfl7l Il8ll . reproduces the exper- 
imentally extracted source function profiles S(r x ), S(r y ) and 
S(r z ). The calculations also show that the source function ex- 
cluding the contribution of resonances (open triangles) is nar- 
rower than the experimentally observed Gaussian. However, 
they do not allow us to draw a firm conclusion concerning the 
value of parameter a. Besides the Therminator default a=Q, 
we tested the value a=— 0.5, the same as used in |@] to describe 
the pion data. Our simulations with a=— 0.5 and the other 
parameters fixed, underestimate the source function S(r z ) al- 
ready for r>5 fm but do not show any change in S(r x ) and 
S(r y ). Substantial improvement can be achieved if we allow 
at the same time Xo to increase from 8 fm/c to ~ 10.5 fm/c. The 
latter value is however considerably bigger than Xo=8.5 fm/c 
reported in [9] for the pions. Given these uncertainties, the 
scenario when kaon freeze-out occurs in the source element 
rest frame from a hypersurface devoid of any space-time cor- 
relation (a=0) is only marginally favoured over the one where 
the emission occurs from the outer surface of the fireball in- 
wards (a<0). 



6 




i i I i i i I i i i i 

f STAR kaon (0-30%) 

STAR kaon (0-20%) 

1 PHENIX kaon (0-30%) 




i i i I i i i i 
HKM Glauber (0-30%) 
Buda-Lund (0-30%) 




Year 2004+2007 
Centrality 0-20% 
k T [GeV/c] 0.2-0.36 



2004 

0-30% 

0.2-0.36 0.36-0.48 



0.6 0.7 0.8 0.9„ 1 

m T (GeV/c 2 ) 



1.1 1.2 



FIG. 4. Transverse mass dependence of Gaussian radii (a) R ou t, (b) 
R s ide an d ( c ) Rlong f° r mid-rapidity kaon pairs from the 30% most 
central Au+Au collisions at ,75/^=200 GeV. STAR data are shown 
as solid stars; PHENIX data llfJl as solid circles (error bars include 
both statistical and systematic uncertainties). Hydro-kinetic model 
(23ll with initial Glauber condition and Buda-Lund model I22I1 cal- 
culations are shown by solid squares and solid curves, respectively. 
The dotted line between the solid squares is to guide the eye. For 
comparison purposes, we also plot the result from the 20% most cen- 
tral Au+Au collisions as open stars. Panel (d) shows corresponding 
experimental values of the Gaussian fit parameter X. 



Although most of the extracted parameters of the expand- 
ing fireball are consistent with those obtained from two-pion 
interferometry ijH, the 3-D source function shapes for kaons 
and pions are very different. This is illustrated in Figure 0(a)- 
(c) which compares the correlation profiles for mid-rapidity 
kaons (circles) with those for mid-rapidity pions (squares) re- 
ported by the PHENIX Collaboration H for the same event 
centrality and transverse momentum selection. The kaon 
source function profiles are generally narrower in width than 
those for pions. Moreover, in contrast to the case for pions, 
a long tail is not observed in the kaon S(r x ) (i.e. along the 
pair's total transverse momentum). Compared to the pion case 
where a prominent cloud of resonance decay pions determines 
the source-function tail profiles in out and long directions @], 
the narrower shape observed for the kaons indicates a much 
smaller role of long-lived resonance decays and/or of the ex- 
ponential emission duration width At on kaon emission. 

Further insight into expanding fireball dynamics can be ob- 
tained by studying the A/r—dependence of the kaon Gaussian 
radii in LCMS. To achieve this goal, in addition to the low- 
est momentum bin (0.2<&t<0.36 GeV/c), we have also anal- 
ysed the kaon pairs with 0.36<&r<0.48 GeV/c. The analysis 



R,[fm] 4.8±0.1±0.2 
R y [fm] 4.3±0.1±0.1 
R z [fin] 4.7±0.1±0.2 

X 0.49±0.02±0.05 
X 2 /ndf 497/289 



4.3±0.1±0.4 
4.0±0.1±0.3 
4.3±0.2±0.4 
0.39±0.01±0.09 
316/283 



4.5±0.2±0.3 
3.7±0.1±0.1 
3.6±0.2±0.3 
0.27±0.01±0.04 
367/283 



TABLE I. Parameters obtained from the 3-D Gaussian source func- 
tion fits for the different datasets. The first errors are statistical, the 
second errors are systematic. 



was carried out for the 30% most central Au+Au collisions at 
\Aaw=200 GeV. This wider centrality cut enabled us to com- 
pare our results to the PHENIX kaon data points obtained at 
higher kj but at the same centrality [10]. A 4-parameter fit to 
the two sets of independent moments with a Gaussian func- 
tion Eq. (|4]i yields a % 2 /ndf of 1.1 and 1.3 respectively. The 
three Gaussian radii and the amplitude obtained from this fit 
are listed in Table Q] Note that the overall normalization of 
S G (r x ,r y ,r : ) may also be affected by systematic factors not 
included in this fit. While the value of X for the 0-20% cen- 
trality data is only marginally smaller than that of Ref. ifioll . 
the analysis of the 30% most central collisions restricted to 
year 2004 data uses looser purity cuts, thus yielding substan- 
tially smaller X. Additional dilution of the correlation strength 
is expected from the (j) — > K + K~ decays, which is, however, 
limited by the (j) decay length of ~1 1 fm in PCMS OH]. Cal- 
culations based on the core-halo model ll20ll employing the 
STAR §/K~ ratio ll2lll yields a maximum 15-20% decrease 
in X at low transverse momenta. Neither of those two effects 
has a significant impact on the values of the extracted Gaus- 
sian radii. 

Figure |4] shows the dependence of the Gaussian radii in 

LCMS (Rout=Rx/Y, Rside=R y and Riong=Rzi Y is tne kine- 
matic Lorentz boost in the outward direction from the LCMS 
to the PCMS frame) as a function of transverse mass wt = 
(m 2 +k\) l l 2 obtained from the fits to the 3-D correlation 
functions from STAR data (stars). The error bars on the 
STAR data are dominated by systematic uncertainties from 
particle identification and momentum resolution. The Gaus- 
sian radii for PHENIX kaon data ifioll (solid circles) are also 
shown, with the error bars representing statistical and system- 
atic uncertainties combined. The model calculations from the 
Buda-Lund model lf22ll and from the Hydro-Kinetic Model 
(HKM) l23tl are shown as solid curves and solid squares, re- 
spectively. While the HKM provides a full microscopic trans- 
port simulation of hydrodynamic expansion of the system fol- 
lowed by dynamic decoupling, the Buda-Lund model is a 
pure analytical solution of the perfect fluid hydrodynamics. 
The latter describes the Gaussian radii of charged pions from 
Au+Au collisions l24ll at the same energy and centrality as 
our kaon data over the whole 0.30<tot<1.15 GeV/c 2 inter- 
val l22tl . Since the exact mx-scaling is an inherent feature of 
perfect fluid hydrodynamics, the Buda-Lund model predicts 
that the kaon and pion radii fall on the same curve. 
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From Figure |4]it is seen that the Gaussian radii for the kaon 
source function display a monotonic decrease with increas- 
ing transverse mass mj from the STAR data at low mj to the 
PHENIX data at higher inj, as do the model calculations of 
Buda-Lund and HKM. The Gaussian radii in the outward and 
sideward directions are adequately described by both models 
over the whole interval. However, there is a marked differ- 
ence between the HKM and the Buda-Lund predictions in the 
longitudinal direction, with the deviation becoming promi- 
nent for mx<0.7 GeV/c 2 where the new STAR data reside. 
Our measurement at 0.2<kj <0.36 GeV/c clearly favours the 
HKM model as more representative of the expansion dynam- 
ics of the fireball, despite the fact that the Buda-Lund model 
describes pion data in all three directions. Hence, exact mj- 
scaling of the Gaussian radii in the longitudinal direction be- 
tween kaons and pions is not supported by our measurements. 

In summary the STAR Collaboration has extracted the 3-D 
source function for mid-rapidity, low transverse momentum 
kaon pairs from central Au+Au collisions at iJsnn=200 GeV 
via the method of Cartesian surface-spherical harmonic de- 
composition. The source function is essentially a 3-D Gaus- 
sian in shape. Comparison with Therminator model calcula- 
tions indicates that kaons are emitted from a fireball whose 
transverse dimension and lifetime are consistent with those 



extracted with two-pion interferometry. However, the 3-D 
source function shapes for kaons and pions are very different. 
The narrower shape observed for the kaons indicates a much 
smaller role of long-lived resonance decays and/or of the ex- 
ponential emission duration width Al on kaon emission. The 
Gaussian radii for the kaon source function display a mono- 
tonic decrease with increasing transverse mass mj over the 
interval 0.55<'«t<1-15 GeV/c 2 . In the outward and side- 
ward directions, this decrease is adequately described by mj- 
scaling. However, in the longitudinal direction, the scaling is 
broken. The results are in favor of the hydro-kinetic predic- 
tions l23tl over pure hydrodynamical model calculations. 
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